/*
 [auto_generated]
 libs/numeric/odeint/test/resize.cpp

 [begin_description]
 This file tests the resize function of odeint.
 [end_description]

 Copyright 2009-2012 Karsten Ahnert
 Copyright 2009-2012 Mario Mulansky

 Distributed under the Boost Software License, Version 1.0.
 (See accompanying file LICENSE_1_0.txt or
 copy at http://www.boost.org/LICENSE_1_0.txt)
 */

// disable checked iterator warning for msvc
#include <boost/config.hpp>
#ifdef BOOST_MSVC
#pragma warning(disable : 4996)
#endif

#define BOOST_TEST_MODULE odeint_resize

#include <vector>
#include <cmath>

#include <boost/test/unit_test.hpp>

#include <boost/numeric/odeint/util/resize.hpp>
#include <boost/fusion/include/vector.hpp>
#include <boost/fusion/include/at.hpp>
#include <boost/units/systems/si/length.hpp>
#include <boost/units/systems/si/time.hpp>
#include <boost/units/systems/si/velocity.hpp>
#include <boost/units/systems/si/acceleration.hpp>
#include <boost/units/systems/si/io.hpp>

using namespace boost::unit_test;
using namespace boost::numeric::odeint;

BOOST_AUTO_TEST_CASE(test_vector) {
  std::vector<double> x1(10);
  std::vector<double> x2;
  resize(x2, x1);
  BOOST_CHECK(x2.size() == 10);
}

BOOST_AUTO_TEST_CASE(test_fusion_vector_of_vector) {
  typedef boost::fusion::vector<std::vector<double>, std::vector<double> > state_type;
  state_type x1;
  boost::fusion::at_c<0>(x1).resize(10);
  boost::fusion::at_c<1>(x1).resize(10);
  state_type x2;
  resize(x2, x1);
  BOOST_CHECK(boost::fusion::at_c<0>(x2).size() == 10);
  BOOST_CHECK(boost::fusion::at_c<1>(x2).size() == 10);
}

BOOST_AUTO_TEST_CASE(test_fusion_vector_mixed1) {
  typedef boost::fusion::vector<std::vector<double>, double> state_type;
  state_type x1;
  boost::fusion::at_c<0>(x1).resize(10);
  state_type x2;
  resize(x2, x1);
  BOOST_CHECK(boost::fusion::at_c<0>(x2).size() == 10);
}

BOOST_AUTO_TEST_CASE(test_fusion_vector_mixed2) {
  typedef boost::fusion::vector<double, std::vector<double>, double> state_type;
  state_type x1;
  boost::fusion::at_c<1>(x1).resize(10);
  state_type x2;
  resize(x2, x1);
  BOOST_CHECK(boost::fusion::at_c<1>(x2).size() == 10);
}

BOOST_AUTO_TEST_CASE(test_fusion_quantity_sequence) {
  namespace units = boost::units;
  namespace si = boost::units::si;

  typedef double value_type;
  typedef units::quantity<si::time, value_type> time_type;
  typedef units::quantity<si::length, value_type> length_type;
  typedef units::quantity<si::velocity, value_type> velocity_type;
  typedef boost::fusion::vector<length_type, velocity_type> state_type;

  state_type x1, x2;
  resize(x2, x1);
}
